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Abstract 


The  report  summarizes  a  research  program  focused  on  the  understanding  of  the  funda¬ 
mental  physics  and  practical  applications  of  Electromagnetic  Local  Flow  Control  (ELEC)  in 
high  speed  flows.  ELEC  encompasses  a  wide  variety  of  energy  deposition  techniques  {e.g., 
plasma  arcs,  laser  pulse,  microwave,  electron  beam,  glow  discharge,  etc)  both  with  and  with¬ 
out  magnetic  held,  and  a  wide  range  of  applications  {e.g.,  drag  reduction,  lift  and  moment 
enhancement,  improved  combustion  and  mixing,  modihcation  of  shock  structure,  etc.).  This 
report  summarizes  two  main  research  achievements  in  the  area  of  ELEC.  The  hrst  is  a  study 
of  the  efficiency  and  effectiveness  of  energy  deposition  in  reducing  drag  on  a  blunt  body, 
and  the  concommitant  effect  on  heat  transfer  to  the  body.  A  computational  study  of  the 
interaction  of  energy  deposition  on  a  blunt  cylinder  was  performed  for  varying  hlament  pulse 
periods.  The  results  indicate  that  energy  deposition  is  both  effective  in  reducing  drag  on 
blunt  bodies  in  supersonic  how,  and  also  efficient  in  the  sense  that  the  power  required  is 
signihcantly  less  than  the  power  saved.  Moreover,  no  adverse  ehects  on  heat  transfer  to 
a  blunt  body  were  observed,  and  in  one  case  a  signihcant  reduction  in  heat  transfer  was 
achieved.  The  second  is  the  development  of  a  fully  three-dimensional,  time-dependent  gas 
dynamic  code  for  microwave  energy  deposition  in  air  incorporating  detailed  kinetics  and 
thermochemistry.  This  code  is  validated  by  means  of  a  comparison  to  experimental  data, 
and  is  used  to  analyze  interaction  of  a  plasma  hlament  with  a  hemisphere  cylinder  and  a 
hemisphere-cone  cylinder.  To  the  authors’  knowledge,  this  is  the  only  code  with  the  capa¬ 
bility  for  simulating  both  the  full  thermochemistry  of  microwave  discharge  in  air  and  the 
interaction  of  the  microwave-generated  plasma  with  an  aerodynamic  shape. 
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I.  Introduction 


In  recent  years  there  has  been  intense  research  activity  in  developing  fundamental  un¬ 
derstanding  and  practical  applications  of  Electromagnetic  Local  Flow  Control  (ELEC).  The 
term  Electromagnetic  Local  Flow  Control  encompasses  a  wide  range  of  energy  deposition 
techniques  including  beamed  energy  deposition  (he.,  electron  beam,  laser,  and  microwave, 
in  the  sense  that  each  represents  beamed  electromagnetic  radiation)  and  DC  discharge  (in 
the  sense  of  energy  deposition  by  means  of  an  imposed  electric  held),  either  with  or  without 
external  magnetic  helds.  This  interest  is  rehected  in  numerous  conferences  and  workshops 
including  the  Weakly  Ionized  Gas  Workshops  (since  1997),  the  St.  Petersburg  Workshops 
(since  2000)  and  the  IVTAN  Workshops  (since  1999).  ELEC  research  encompasses  a  wide 
variety  of  energy  deposition  techniques  {e.g.,  plasma  arcs,  laser  pulse,  microwave,  electron 
beam,  glow  discharge,  etc)  both  with  and  without  magnetic  held,  and  a  wide  range  of  ap¬ 
plications  {e.g.,  drag  reduction,  lift  and  moment  enhancement,  improved  combustion  and 
mixing,  modihcation  of  shock  structure,  etc.).  Reviews  of  this  research  include  Shang,^ 
Knight  et  al,'^  Fomin  et  al,^  Knight,'^  Bletzinger  et  al^  and  Knight.® 

The  advantages  of  Electromagnetic  Local  Flow  Control  are  twofold.  First,  the  actua¬ 
tion  times  are  comparable  or  substantially  smaller  than  the  howheld  development  times  in 
supersonic  and  hypersonic  hight.  For  example,  consider  a  hypersonic  vehicle  operating  at 
Mach  6  and  30  km  altitude.  The  relevant  timescale  for  howheld  development  in  the  vicinity 
of  an  /  =  10  cm  dimension  {e.g.,  a  cowl  leading  edge  radius)  is  //Goo  =  5.5  ■  10“®  s.  Typical 
pulse  durations  for  laser  and  microwave  are  0(10“®  s)  and  0(10“®  s),  respectively,  thereby 
permitting  modihcation  of  the  local  howheld.  Electromechanical  or  mechanical  how  control 
systems,  however,  are  much  slower  with  characteristics  actuation  times  that  are  0(10“^  s). 
Second,  beamed  energy  deposition  can  be  ehectively  used  to  modify  the  howheld  at  large 
distances  from  the  energy  source  {i.e.,  at  distances  from  the  vehicle),  whereas  electromechan¬ 
ical  or  mechanical  how  control  systems  are  limited  to  modifying  the  how  in  the  immediate 
vicinity  of  the  vehicle. 

Plasmadynamics  is  an  inherent  element  in  Electromagnetic  Local  Flow  Control.  Typi¬ 
cal  methods  of  energy  deposition  in  air  by  electromagnetic  techniques  {e.g.,  DC  discharge, 
electron  beam,  laser  and  microwave)  generate  a  high  temperature  plasma  with  ions  (N^, 
O/,  NO’*',  N+,  etc),  electrons  (e),  dissociated  species  (N,  O),  excited  states  {N2{A^T,u+), 
N2{B^Ilg) ,  etc)  and  molecules  (NO)  in  addition  to  N2  and  O2.  The  formation  of  the  plasma 
{e.g.,  its  shape  and  structure)  has  an  important  ehect  on  its  ability  to  modify  the  howheld. 

For  example,  microwave  energy  deposition  in  air  gener¬ 
ates  a  plasma  characterized  by  a  “warm”  region  whose 
size  is  on  the  order  of  the  microwave  wavelength  and 
which  contains  thin  “hot”  hlaments  whose  diameter 
is  one  to  two  orders  of  magnitude  smaller  (Fig.  1). 

Detailed  experiments  and  computations  have  demon¬ 
strated  that  the  interaction  of  the  microwave-generated 
plasma  with  a  blunt  body  in  supersonic  how  can  lead  to 
a  momentary  reduction  in  drag  principally  due  to  the 
interaction  of  the  hlaments  with  the  blunt  body  shock 

(Kolesnichenko  et  and  Wilkerson  et  aA^).  Figure  1.  Microwave  filaments 
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This  report  summarizes  two  main  research  achievements  in  the  area  of  ELFC.  The  first 
is  a  study  of  the  efficiency  and  effectiveness  of  energy  deposition  in  reducing  drag  on  a  blunt 
body,  and  the  concommitant  effect  on  heat  transfer  to  the  body.  A  computational  study  of 
the  interaction  of  energy  deposition  and  a  blunt  cylinder  is  performed  for  varying  filament 
pulse  periods.  This  study  is  detailed  in  Section  II.  The  second  is  the  development  of  a  fully 
three-dimensional,  time-dependent  gas  dynamic  code  for  microwave  energy  deposition  in  air 
incorporating  detailed  kinetics  and  thermochemistry.  This  code  is  validated  by  means  of  a 
comparison  to  experimental  data,  and  is  used  to  analyze  interaction  of  a  plasma  hlament 
with  a  hemisphere  cylinder  and  a  hemisphere-cone  cylinder.  These  results  are  presented  in 
Section  III. 

II.  Effectiveness,  Efficiency  and  Heat  Transfer  Effects  of 
Energy  Deposition  for  Flow  Control 

A  computational  study  was  performed  to  examine  the  interaction  of  a  periodic  sequence 
of  heated  hlaments  with  a  blunt  cylinder  in  supersonic  flow.  The  heated  hlaments  are 
introduced  on  the  axis  of  symmetry  upstream  of  the  cylinder  and  aligned  with  the  freestream 
flow.  The  heated  hlaments  are  an  idealization  of  the  hlamentary  plasma  structures  generated 
by  a  microwave  discharge  in  air  (Fig.  1).  The  interaction  is  examined  with  a  focus  on 
two  effects,  namely,  1)  the  efficiency  and  effectiveness  of  drag  reduction,  and  2)  the  effect 
on  heat  transfer  to  the  cylinder  face.  The  conhguration  is  shown  in  Fig.  2.  The  value 
of  the  freestream  Reynolds  number  Re  =  PooUooD j Hoo  is  chosen  to  approximate  typical 
experimental  values,  such  as  experiments  performed  by  Kolesnichenko  et  al^‘^  where  typical 
values  for  poo,  Tx)  and  D  are  25-50  Torr,  155  K,  and  20-40  mm,  respectively.  The  how 
parameters  are  dehned  in  Table  1.  Note  that  the  hxed  cylinder  surface  temperature 
exceeds  the  adiabatic  wall  temperature  Ta  by  11%  (hot  wall).  The  pulse  period  L/D  is 
varied  from  4/3  to  inhnitely  long. 


Figure  2.  Computational  domain 


The  heated  hlaments  are  dehned  by  the  density  ratio  a  =  Pf/poo  (where  pf  is  the  hlament 
density),  diameter  ratio  d/ D,  length  ratio  1/ D  and  period  Lj D.  Fixed  values  of  a  =  0.5  and 
l/D  =  1  (except  for  the  inhnitely  long  hlament)  were  assumed.  The  hlaments  are  introduced 
at  the  inhow  boundary  of  the  computational  domain  with  static  pressure  pf  =  Poo  and 
velocity  Uf  =  I/qo-  The  temperature  of  the  hlament  is  therefore  Tj  =  a~^T^. 

The  governing  equations  are  the  unsteady  laminar  compressible  Navier  Stokes  equations 
for  an  ideal  gas.  The  how  is  axisymmetric.  The  molecular  viscosity  is  dehned  by  Sutherland’s 
Law.  A  code  was  written  in  C-I--I-  to  solve  the  governing  equations  using  a  hnite  volume 


4  of  17 


Table  1.  Flow  Parameters 


Type 

Description 

Definition 

Value 

Dimensionless  Flow 

Mach 

Moo 

1.89 

Specific  heat  ratio 

7 

1.4 

Reynolds  number 

Re 

o 

o 

Prandtl  number 

Pr 

0.7 

Dimensionless  Filament 

Density  ratio 

a 

0.5  all  cases 

Diameter 

d 

D 

0.1  all  cases 

Length 

1 

D 

1.0,  1.0,  1.0,  cx) 

Pulse  (Period) 

L 

D 

|,  2.0,  4.0,  CX) 

Cylinder 

Temperature 

Tw/Too 

2.0  all  cases 

Adiabatic  Wall  Temperature 

Ta/Too 

1.8 

method.  Details  of  the  code  are  given  in  Anderson.  Van  Leer’s  flux  splitting  method^^ 
is  used  for  the  inviscid  fluxes.  The  viscous  fluxes  are  solved  with  a  second  order  accurate 
discretization.  The  conservative  variables  are  reconstructed  using  a  third  order  accurate 
Essentially  Non-Oscillatory  (ENO)  Method. A  second  order  accurate  explicit  Runge-Kutta 
scheme  is  used  to  integrate  in  time.^®  The  code  is  parallelized  using  Message  Passing  Interface 
(MPI).  The  corner  of  the  cylinder  creates  a  local  singularity  to  the  flow.  To  minimize  errors 
introduced  by  the  singularity,  an  entropy  £x  proposed  by  Woodward  et  al^'^  was  incorporated 
to  the  code.  For  the  larger  period  cases  {L/D  =  2  and  L/D  =  4),  the  density  and  pressure 
were  relaxed  in  time  if  the  computed  density  in  the  vicinity  of  the  corner  was  momentarily 
negative.  The  density  and  pressure  were  set  equal  to  a  weighting  between  the  computed 
parameters  from  the  current  timestep  to  the  computed  parameters  from  the  previous  step 
in  the  Runge-Kutta  scheme. 

The  computational  domain  is  shown  in  Fig.  3.  The  pulses  have  the  greatest  effect  at  the 
face  of  the  cylinder.  Therefore,  only  one  cylinder  diameter  in  length  of  the  body  was  included 
in  the  flow  domain.  The  height  of  the  domain  extends  two  and  a  half  cylinder  diameters 
above  the  body.  The  blunt  body  shock  is  expected  to  lense  forward  upon  interaction  with 
the  heated  hlaments  based  on  trial  runs,  and  consequently  the  upstream  boundary  of  the 
flow  domain  is  two  cylinder  diameters  in  length  in  front  of  the  body. 

The  principal  measures  of  aerodynamic  performance  associated  with  the  interaction  of 
the  heated  filaments  with  the  blunt  cylinder  are  the  efficiency  and  ejfectiveness  of  the  drag 
reduction.  The  efficiency  is  dehned  as  the  ratio  of  the  power  saved  due  to  drag  reduction 
by  the  heated  filament  to  the  power  necessary  to  create  the  heated  hlament.  The  efficiency 
is  obtained  by  dividing  the  average  thrust  power  saved  due  to  the  hlament  divided  by  the 
power  used  to  create  the  filament.  The  efficiency  of  drag  reduction  in  terms  of  dimensionless 
variables  is 


V 


8(7  -  mi 

1  —  a 


D  DA 
/  Jo  Jo 


{Po  —  p)rdrdt 


(1) 


where  the  pressure  is  nondimensionalized  by  pooU^,  the  subscript  o  indicates  the  absence 
of  the  filament  and  l/D  =  L/D  =  Too  for  L/D  =  oo.  Derivation  of  this  expression  may  be 
found  in  Anderson. 

The  effectiveness  is  dehned  as  the  ratio  of  the  average  frontal  drag  reduction  due  to  the 
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Figure  3.  Axisymmetric  wedge  solution  domain 


presence  of  the  filament  to  the  frontal  drag  in  the  absence  of  the  hlament.  The  effectiveness 

is 

.  ^  fo^^jPo  -  p)rdrdt 
L/D  Pordr 

where  L/D  =  t^o  for  L/D  =  oo. 

The  principal  measnre  of  aerothermodynamic  performance  associated  with  the  interaction 
of  the  heated  hlaments  with  the  blunt  cylinder  is  the  ratio  <h  of  the  average  heat  transfer 
on  the  cylinder  face  due  to  the  interaction  with  the  hlaments  to  the  heat  transfer  on  the 
cylinder  face  in  the  absence  of  the  hlaments. 


$  = 

Qnofil 


(3) 


where 


_  1  r  /•1/2 

Q=  /  q^rdrdt 


(4) 


r  JO  JO 

where  r  is  the  time  interval  non-dimensionalized  by  D/Uoo  at  which  the  integral  is  taken, 
r  is  the  radius,  non-dimensionalized  by  the  cylinder  diameter,  and  is  the  heat  hux  to 
the  cylinder  face  non-dimensionalized  by  pooU^.  A  positive  heat  hux  implies  cooling  of  the 
cylinder  face.  The  sensitivity  of  the  heat  transfer  ratio  to  the  grid  size  was  determined. 
For  each  pulse  period  evaluated,  two  grids  were  allowed  to  reach  a  statistically  stationary 
state,  the  second  with  cell  sizes  half  the  size  of  the  initial.  The  error  of  the  rehned  grid  is 
calculated  based  on  these  two  cases  to  determine  the  accuracy  of  the  computations.  The 
estimated  uncertainty  in  $  in  the  rehned  grid  is  less  than  6%. 

The  efficiency  rj,  ehectiveness  (  and  heat  transfer  parameter  <F  vs  pulse  period  L/D  are 
presented  in  Figs.  4,  5,  and  6,  respectively.  The  overall  efficiency  is  very  high  with  typical 
values  approximately  100  (he.,  for  every  Watt  of  power  supplied  to  heat  the  hlament,  more 
than  100  Watts  of  thrust  power  are  saved)  with  a  maximum  value  at  L/D  =  2.  The 
extraordinarily  large  values  of  efficiency  are  consistent  with  the  analysis  of  Kolesnichenko  et 
al}^  Substantial  levels  of  drag  reduction  up  to  30%  are  achieved  with  pulsed  hlaments.  The 
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effectiveness  decreases  with  increasing  period  as  expected  since  fewer  filaments  interact  with 
the  blnnt  cylinder  per  nnit  time. 

A  major  concern  is  the  effect  of  the  heated  filaments  on  the  heat  transfer  to  the  cylinder 
face.  It  is  notable  that  the  heat  transfer  parameter  $  does  not  exceed  nnity  implying  that 
there  is  no  increase  in  average  heat  transfer  to  the  blnnt  cylinder  face.  Indeed,  the  heat 
transfer  is  reduced  for  the  L/D  =  4/3  case.  The  reason  for  the  minimnm  heat  transfer  at 
L/D  =  4/3  is  discnssed  below  in  terms  of  the  vortex  strnctnre. 


Figure  4.  Efficiency  Figure  5.  Effectiveness  Figure  6.  Heat  Transfer 


Contours  of  velocity  magnitude  and  instantaneous  streamlines  are  shown  in  Fig.  7  for 
L/D  =  oo.  The  interaction  of  the  filament  with  the  blunt  body  shock  generates  an  un¬ 
stable  shear  layer  (identihed  by  the  abrupt  change  in  color  from  blue  [low  speed]  to  green 
[high  speed]  in  Fig.  7)  resulting  in  the  formation  of  sequence  of  periodic  vortices.  As  the 
line  of  vortices  interacts  with  the  cylinder  face,  two  quasi-stationary  vortical  structures  are 
established.  A  counter-clockwise  rotating  vortex  forms  just  below  the  line  of  vortices  (CCW 
vortex)  and  a  clockwise  rotating  vortex  (CW  vortex)  develops  beneath  it  (see  Farzan  et  al}^ 
The  size  and  location  of  the  vortices  vary  in  time,  but  their  presence  ahead  of  the  cylinder 
face  is  maintained.  Instantaneous  streamlines  from  the  freestream  above  the  high  temper¬ 
ature  filament  remain  above  the  line  of  vortices.  This  deflection  of  streamlines  creates  a 
perceived  streamlining  of  the  body  and  is  responsible  for  the  reduction  in  drag. 

The  pulsed  hlament  case  with  a  period  of  L/D  =  4/3  demonstrates  the  greatest  reduction 
in  average  heat  transfer  to  the  cylinder  face  (Fig.  6).  The  contours  of  the  velocity  magnitude 
and  the  instantaneous  streamlines  are  shown  in  Fig.  8.  Each  pulsed  hlament  generates  a  vor¬ 
tex  behind  the  bow  shock.  The  proximity  of  the  vortices  is  sufficient  to  provide  streamlining 
for  the  majority  of  the  heat  transfer  cycle,  which  creates  a  stagnation  region  below  the  line 
of  vortices.  The  nearly  stagnant  how  warms  to  the  cylinder  temperature,  creating  a  thermal 
layer  at  the  cylinder  face.  This  quasi-stationary  thermal  layer  then  decreases  heat  transfer 
from  the  cylinder  surface.  The  vortices  established  in  the  inhnitely  long  hlament  case  are 
not  generated  in  this  case.  This  is  likely  due  to  a  brief  break  between  streamlining  of  the 
vortices,  which  sweeps  how  ahead  of  the  cylinder  over  the  corner.  Thus,  the  heat  transfer  is 
reduced  further  for  this  case  than  the  inhnitely  long  hlament  case. 

The  velocity  magnitude  contours  and  instantaneous  streamlines  for  the  periodic  energy 
deposition  case  with  a  period  of  L/D  =  2  are  shown  in  Fig.  9.  The  integrated  heat  transfer 
is  close  to  unity  for  this  case.  The  vortices  provide  a  partial  streamlining,  which  allows  for 
a  heated  layer  to  build  ahead  of  the  cylinder  face.  The  time  between  hlaments,  however,  is 
sufficiently  large  that  the  heated  layer  is  convected  away  from  the  body,  thereby  introducing 
cooler  how  into  the  region  ahead  of  the  cylinder  face  and  increasing  the  heat  transfer.  The 
time  between  streamlining  is  much  greater  than  the  time  the  streamlining  creates  a  thermal 
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Figure  7.  L/D  =  oo 


Figure  8.  LjD  =  4/3 


layer,  thus  the  integrated  heat  transfer  is  not  greatly  reduced. 

The  case  of  the  pulsed  hlament  with  L/D  =  4  exhibits  a  heat  transfer  parameter  <h  1. 
The  velocity  magnitude  contours  and  instantaneous  streamlines  are  shown  in  Fig.  10.  Al¬ 
though  the  vortices  temporarily  streamline  the  body  and  allow  for  partial  warming  of  the 
flow  ahead  of  the  cylinder  face,  the  time  between  vortices  is  large  compared  to  the  time 
wherein  the  vortex  streamlines  the  flow.  Hence,  the  heat  transfer  is  close  to  one,  similar  to 
the  pulsed  period  of  two  case. 


X/D 


X/D 


Figure  9.  L/D  =  2 


Figure  10.  L/D  =  A 
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III.  Interaction  of  Microwave-generated  Plasma  with 
Blunt  Aerodynamic  Shapes 


A  fully  three-dimensional,  time-dependent,  viscous  gas  dynamic  code  for  simulating  mi¬ 
crowave  energy  deposition  in  air  has  been  developed  incorporating  detailed  kinetics  and 
thermochemistry.  The  governing  equations  and  thermochemistry  model  are  described  in 
Knight  et  al.^^  The  code  is  written  in  C-I--I-  and  parallelized  using  Message  Passing  Interface 
(MPI).  To  the  authors’  knowledge,  this  is  the  only  code  with  the  capability  for  simulating  both 
the  full  thermochemistry  of  microwave  discharge  in  air  and  the  interaction  of  the  microwave¬ 
generated  plasma  with  an  aerodynamic  shape.  The  thermochemistry  model  was  developed  in 
the  laboratory  of  Dr.  Yuri  Kolesnichenko  (Joint  Institutes  for  High  Temperatures,  Moscow, 
Russia).  The  code  was  validated  by  comparison  to  experimental  data  for  the  interaction  of 
a  microwave-generated  plasma  with  a  hemisphere-cylinder  and  hemisphere-cone-cylinder. 


A.  Plasma  Interaction  with  Hemisphere  Cylinder 


The  gas  dynamic  code  is  applied  to  the  simulation  of  the  interaction  of  a  microwave-generated 
plasma  with  the  flow  past  a  hemisphere-cylinder  at  Moo  =  2.1.  Fig.  11  displays  the  flow  con- 
hguration  which  is  axisymmetric.  The  microwave  pulse  is  focused  at  a  distance  Xo  upstream 
of  the  center  of  curvature  of  the  hemisphere.  The  freestream  conditions  are  listed  in  Table  2. 
The  inflow  corresponds  to  standard  air  at  the  specihed  static  pressure  and  temperature  with 
the  inclusion  of  a  small  concentration  of  electrons  N^.  and  molecular  oxygen  ions  Nq+.  The 
freestream  concentration  of  all  other  species  is  set  to  zero.  The  microwave  parameters  are 
indicated  in  Table  3. 


Figure  11.  Computational  domain  Figure  12.  Initial  condition  for  MW  pulse 


Details  of  the  grids  are  listed  in  Table  4  and  discussed  in  Knight  et  al.'^^  The  number 
of  cells  and  in  the  f-  and  //-directions,  respectively,  are  indicated,  together  with  the 
total  number  of  cells  iVt^tai  in  the  computational  domain.  The  number  of  processors  Np  =  36. 
Axis  boundary  conditions  are  applied  on  the  symmetry  axis  AB.  Freestream  conditions  are 
imposed  on  the  outer  boundary  BC,  and  zero  gradient  conditions  at  the  outflow  boundary 
CD.  Tangential  (slip)  boundary  conditions  are  applied  on  the  hemisphere-cylinder  surface 
AD  since  the  flow  was  assumed  inviscid  for  this  conhguration.  The  initial  condition  for  the 
microwave  pulse  is  the  converged  solution  for  flow  past  the  body  at  the  freestream  conditions 
shown  in  Table  2.  The  initial  computed  static  pressure  contours  are  shown  in  Figs.  12. 
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Table  2.  Preestream  Conditions 


Table  3.  Microwave  Parameters 


Quantity 

Value 

Quantity 

Value 

2.1 

Eo  (kV/cm) 

2.3 

Poo  (Torr) 

26 

A  (cm) 

3.33 

Too  (K) 

154 

/  (GHz) 

9.0 

D  (cm) 

2.0 

To  (/is) 

1.20 

Nm^  (cm-3) 

1.287-101® 

Ti  (ps) 

1.21 

N02  (cm“3) 

3.422  ■  lOi^ 

Xo  (cm) 

-3.5 

No+  (cm-3) 

3.0  ■  10^ 

Ag  (cm“^) 

0 

0 

Table  4.  Details  of  Grids 


Quantity  Grid  No.  1  Grid  No.  2 


276 

550 

90 

180 

-^^total 

24,840 

99,000 

Ar  ■ 

(cm) 

4.00  ■  10-2 

2.00  ■  10-2 

Ar^,, 

(cm) 

4.00  ■  10-2 

2.00  ■  10-2 

^^min 

(cm) 

1.13  ■  10-2 

0.57-10-2 

A5 

max 

(cm) 

5.86  ■  10-2 

2.95  -  10-2 

Np 

36 

36 

The  development  of  the  pressure  p,  temperature  T  and  electron  temperature  Tg  at  the 
center  of  the  discharge  is  shown  in  Fig.  13.  Also  shown  is  the  imposed  electric  held  Eof{t) 
and  reduced  held  E/M  at  the  center.  The  pressure  reaches  a  maximum  value  of  346  Torr 
at  t  =  1.4  /is,  and  the  gas  temperature  reaches  a  maximum  of  1958  K  at  t  =  1.4  /xs.  The 
maximum  reduced  held  is  135  Td.  The  contributions  to  the  gas  heating  (^elastic,  <?jouie,  and 
^'reactions  at  the  Center  of  the  discharge  are  shown  in  Fig.  14.  The  most  signihcant  contribution 
is  (^reactions  which  rcachcs  a  maximum  value  of  1.25  MW/cm^  at  t  =  1.2  /is.  The  corresponding 
maximum  values  for  and  geiastic  are  0.19  MW/cm^  and  0.019  MW/cm^.  Note  that  the 
rotational  relaxation  factor  v  is  on  the  order  of  10“^  for  a  reduced  held  of  100  Td,  thereby 
accounting  for  the  low  level  of  Joule  heating.  The  Joule  heating  term  becomes  zero  at  the 
termination  of  the  pulse  since  the  electric  held  vanishes.  Since  the  collision  frequency  z/g 
is  an  assumed  function  of  the  reduced  held,  the  elastic  heating  also  goes  to  zero  when  the 
electric  held  vanishes.  Overall,  the  principal  contribution  to  the  gas  heating  is  due  to  the 
thermochemical  reactions  (he.,  the  quenching  of  excited  states). 

Details  of  the  kinetics  of  the  plasma  formation  are  illustrated  in  Figs.  15  and  16  which 
display  the  evolution  of  the  concentration  (cm“^)  of  selected  species  at  the  center  of  the 
plasma  for  f  =  0  to  t  =  2  /is.  Fig.  15  indicates  a  modest  dissociation  of  N2  with  a  peak 
value  of  N  equal  to  3.6%  of  the  freestream  N2  concentration.  The  concentration  of  the 
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electronically  excited  states  of  nitrogen  N(^D),  N{^P),  N2{a^T,~),  N2{A^T,u+),  N2{B^Ilg), 
N2{C^IIu)  do  not  exceed  4  ■  10^®  cm“^  {i.e.,  3.1%  of  the  freestream  N2  concentration)  and 
rapidly  decay  (Figs.  15  and  16). 


E  (V/m) 


Figure  13.  p,  T,  Tg  vs  t 


Figure  14.  Heating  terms 


Figure  15.  Plasma  species  Figure  16.  Plasma  species 

The  computed  and  experimental  surface  pressure  vs  time  on  the  centerline  of  the  body 
surface  for  the  hemisphere-cylinder  are  shown  in  Fig.  17.  The  natural  frequency  of  the 
experimental  pressure  transducer  is  approximately  120  kHz  and  therefore  the  computed 
surface  pressure  was  hltered  to  simulate  the  experimental  measurement.  The  filtering  was 
performed  by  multiplying  the  Fourier  coefficients  of  the  computed  surface  pressure  time  series 
by  a  factor  c(/)  where 


f  1  0</</i 

c{f)=\  1  -  (/  -  /l)/(/2  -  /l)  /l</</2  (5) 

i  0  f2<f 

where  /i  =  30  kHz  and  /2  =  60  kHz.  The  time  point  of  initial  pressure  decrease  was 
matched  between  the  computation  and  experiment  in  order  to  account  for  the  uncertainty  in 
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the  location  of  the  microwave  pulse  in  the  experiment.  This  resulted  in  a  +10/is  shift  in  the 
computed  pressure.  The  upstream  undisturbed  centerline  pressure  levels  were  also  aligned  to 
simplify  comparison  of  the  pressure  variation  between  computation  and  experiment,  implying 
a  correction  of  4  Torr  and  7  Torr  for  Grid  Nos.  1  and  2,  respectively.  The  agreement  between 
the  computed  and  experimental  surface  pressure  is  good.  The  hner  grid  solution  agrees 
more  closely  with  the  experiment.  The  maximum  experimental  pressure  decrease  is  30  Torr, 
and  the  maximum  predicted  pressure  decrease  (Grid  No.  2)  is  23.4  Torr.  The  effect  of 
the  hltering  on  the  computed  centerline  pressure  for  Grid  No.  2  is  shown  in  Fig.  18.  The 
hltering  signihcantly  reduces  the  magnitude  of  the  computed  pressure  drop  associated  with 
the  interaction  of  the  microwave-generated  plasma  and  the  hemisphere  cylinder. 


Figure  17.  Centerline  pressure  Figure  18.  Centerline  pressure 


B.  Plasma  Interaction  with  Hemisphere-Cone  Cylinder 

The  gas  dynamic  code  is  applied  to  the  simulation  of  the  interaction  of  a  microwave-generated 
plasma  with  the  flow  past  a  hemisphere-cone-cylinder.  The  flow  conhguration  is  shown  in 
Fig.  19  and  corresponds  to  experiments.^®  The  microwave  pulse  is  focused  at  a  distance  Xo 
upstream  of  the  center  of  curvature  of  the  hemisphere.  The  freestream  conditions  are  shown 
in  Table  5  and  the  microwave  parameters  in  Table  6.  The  freestream  species  concentrations 
for  N2  and  O2  correspond  to  standard  air  at  the  specihed  static  pressure  and  temperature.  A 
small  freestream  concentration  of  electrons  Ne  and  molecular  oxygen  ions  Nq+  were  included. 
The  freestream  concentration  of  all  other  species  is  set  to  zero.  The  Reynolds  number  based 
upon  the  cylinder  diameter  and  freestream  conditions  is  Rcoo  =  2.47  ■  10®.  The  Prandtl 
number  Pr  =  0.72  and  the  Lewis  number  Le  =  1. 

Two  separate  computations  were  performed  (Gase  Nos.  1  and  2)  with  different  shapes  of 
the  electric  held  (Table  6).  The  cases  differed  in  two  respects.  First,  the  streamwise  length  of 
the  electric  held  is  equal  to  |A  and  A  for  Gase  Nos.  1  and  2,  respectively.  This  is  the  principal 
diherence  between  the  two  cases.  Second,  Gase  No.  1  is  inviscid  (since  it  was  a  preliminary 
computation  to  assess  the  ehect  of  If),  and  Gase  No.  2  is  viscous.  The  number  of  cells 
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shock 
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Figure  19.  Computational  domain 


Table  5.  Freestream  Table  6.  Microwave 


Table  7.  Details  of  Computation 


Quantity 

Value 

Quantity 

Case  1 

Case  2 

Moo 

2.1 

Eo  (kV/cm) 

3.35 

3.35 

Poo  (Torr) 

41 

A  (cm) 

3.33 

3.33 

Too  (K) 

153 

/  (GHz) 

9.0 

9.0 

(K) 

288 

To  (/is) 

1.20 

1.20 

D  (cm) 

4.0 

Ti  (/is) 

1.21 

1.21 

Nn2  (cm-3) 

2.044  ■  10^® 

Xo  (cm) 

-2.1 

-2.1 

N02  (cm-3) 

5.434  ■  10^^ 

If  (cm) 

1.67 

3.33 

No+  (cm-3) 

0 

0 

n 

0 

1 

Nf.  (cm-^) 

0 

0 

Model 

I 

V 

Quantity 


Nr, 

Ar  ■ 

(cm) 

Ar^,, 

(cm' 

^^min 

(cm) 

A5 

(cm 

Np 

Value 

350 

286 

100,100 
1.00  ■  10-2 
1.00  ■  10-2 
1.00  ■  10-2 
7.19-10-2 
44 


LEGEND 

I  Inviscid  V  Viscous 


and  in  the  Q  and  ?7-directions,  respectively,  are  indicated,  together  with  the  total  number 
of  cells  iVtotai  in  the  computational  domain  in  Table  7.  The  grid  spacing  Ar  in  the  r^-direction 
is  uniform  and  equal  to  0.0025/1.  The  minimum  grid  spacing  As  in  the  .^-direction  on  the 
hemisphere  cylinder  is  0.0025/1  and  the  maximum  As  is  0.018/1.  The  number  of  processors 
Np  =  44.  Axis  boundary  conditions  are  applied  on  the  symmetry  axis  AB.  Freestream 
conditions  are  imposed  on  the  outer  boundary  BC,  and  zero  gradient  conditions  at  the 
outflow  boundary  CD.  For  Case  1,  slip  boundary  conditions  are  applied  on  the  hemisphere- 
cone-cylinder  surface  AD.  For  Case  2,  no-slip  isothermal  boundary  conditions  are  applied 
on  the  hemisphere-cone-cylinder  surface  AD  with  =  I.88T00.  The  semi-implicit  Runge- 
Kutta  algorithm  is  used  in  the  initial  stages  of  the  computation  due  to  the  rapid  increase 
in  species  concentrations  and  gas  temperature  in  the  region  of  the  microwave  beam  during 
the  microwave  pulse.  After  a  few  microseconds  of  physical  time,  the  temporal  algorithm  is 
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switched  to  an  explicit  Runge-Kutta  scheme.  The  initial  condition  for  the  microwave  pnlse 
is  the  converged  solntion  for  flow  past  the  body  at  the  freestream  conditions. 


Figure  20.  Case  1  Figure  21.  Case  2 

The  computed  and  experimental  surface  pressure  vs  time  on  the  centerline  of  the  body 
surface  for  the  hemisphere-cone-cylinder  are  shown  in  Figs.  20  and  21  for  Case  Nos.  1  and  2, 
respectively.  Due  to  the  experimental  uncertainty  in  the  location  of  the  microwave  discharge, 
the  timescale  of  the  computation  is  shifted  10.89  fvs  and  19.89  fj,s  for  Case  Nos.  1  and  2, 
respectively,  to  align  the  initial  pressure  rise  with  the  experiment.  The  predicted  minimum 
pressure  in  Case  1  shows  excellent  agreement  with  experiment.  However,  it  is  evident  that 
the  duration  of  the  interaction  is  too  short  compared  with  the  experiment.  This  implies 
that  the  streamwise  length  of  the  microwave-generated  plasma  is  too  small.  The  streamwise 
length  of  the  plasma  is  sensitive  to  the  streamwise  shape  (x— direction  shape)  of  the  electric 
held.  For  Case  1,  If  =  ^A,  and  therefore,  a  greater  streamwise  length  If  =  \  of  the  electric 
held  was  assumed  for  Case  2.  The  predicted  minimum  pressure  in  Case  2  shows  excellent 
agreement  with  experiment.  Moreover,  the  duration  of  the  pressure  drop  more  closely  agrees 
with  experiment.  Nonetheless,  further  research  is  needed  to  dehne  a  more  accurate  spatial 
distribution  of  the  electric  held. 

Examples  of  the  instantaneous  how  structure  are  presented  in  Figs.  22  and  23.  Temper¬ 
ature  contours  are  displayed  in  Fig.  22  at  t  =  28.7  fj,s  corresponding  to  the  peak  centerline 
surface  pressure  in  Fig.  21.  The  time  indicated  in  the  hgures  includes  the  ohset  of  19.89  /xs 
introduced  in  Fig.  21  to  align  the  experimental  and  computed  peak  pressure  associated  with 
the  impingement  of  the  blast  wave  on  the  centerline  surface  of  the  hemisphere.  The  inter¬ 
action  of  the  heated  plasma  causes  a  lensing  upstream  of  the  blunt  body  shock.  Contours 
of  concentration  of  atomic  oxygen  (formed  by  dissociation  of  O2  during  the  microwave  dis¬ 
charge)  are  shown  in  Fig.  23  together  with  contours  of  static  pressure  at  the  same  instant  as 
Fig.  22.  High  levels  of  atomic  oxygen  are  present,  thus  indicating  the  importance  of  modeling 
the  flowfield  chemistry  associated  with  the  microwave  discharge. 
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Figure  22.  T  at  t  =  28.7  //s 


Figure  23.  Nq  and  p  at  t  =  28.7  fis 
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